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The two-fluid and Landau criteria for superfluidity are compared for trapped Bose gases. While 
the two-fluid criterion predicts translational superfluidity, it is suggested, on the basis of the ho- 
mogeneous Gross-Pitaevski limit, that a necessary part of Landau's criterion, adequate for non- 
translationally invariant systems, does not hold for trapped Bose gases in the GP limit. As a 
consequence, if the compressibility is detected to be very large (infinite by experimental standards), 
the two-fluid criterion is seen to be the relevant one in case the system is a translational superfluid, 
while the Landau criterion is the relevant one if translational superfluidity is absent. 
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The experimental observation of Bose Einstein conden- 
sation (BEC) in a vapor of Rb atoms at very low tem- 
peratures (of the order of 170 x 10~^K) in 1995 [1] chal- 
lenged theoretical and mathematical physicicists to find 
an explanation of the phenomenon. Since atom-atom in- 
teractions are weak - the range of interaction was about 
10~^cm, while at the required densities the interatomic 
spacing was about 10~''cm- the gas could be considered 
to be almost perfect except for elastic collisions of hard 
core type, and a theory was soon developed in terms of 
the so-called Gross Pitaevski (G.P.) functional or equa- 
tion [2] which may be derived in the G.P. limit, defined 
by the condition: 
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where a denotes the scattering length, N is the number of 
particles and V = is the volume enclosing the system. 
A mathematical proof of BEC was only achieved in 2002 
in a remarkable paper by Lieb and Seiringer. [3]. 

Setting /i = /2m , the Hamiltonian describing the 
system may be written 
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where d is a positive potential satisfying certain regular- 
ity conditions, and Hq is the one-particle Hamiltonian 
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acting on L^(R'^, daf ) , with A the Laplacian on R'^. 
Above, Hq describes particles confined in a trap poten- 
tial Vt (such that Vt{x) ^ oo as ^ oo) (see, 
e.g., ref. [4]). In this letter, we address the (not yet en- 
tirely clear) issue of superfluidity, most particularly that 
of translational superfluidity. In [5] the following crite- 
rion was used to prove superfluidity of these gases. Let 
Eq denote the ground state energy of the system in the 



rest frame and Eq the ground state energy, measured in 
the moving frame, when a velocity field v is imposed. 
Then for small v 
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where N is the particle number and m the mass. The er- 
ror term must be bounded independently of N (see [12], 
pg.44) and ps and p denote, respectively, the density of 
the superfluid and the total particle density in the two- 
fluid picture (we shall call this "the two-fluid criterion" 
for brevity). By (4) the free Bose gas in the ground state 
is a superfluid, with Eq — 0, ps — p and Eq — \/2Nm{v)'^ 
and no error term, while it is not one by the Landau crite- 
rion, and thus the two criteria differ in principle. The free 
Bose gas is not a physical system, but now that trapped 
Bose gases are an experimental reality, the fundamen- 
tal question poses itself a) do both the two-fluid and the 
Landau criterion predict translational superfluidity for 
these systems and b) if not, which criterion describes the 
physical property of superfluidity correctly? In order to 
address this question we must now turn to a brief analysis 
of the Landau criterion. 

The original Landau criterion was formulated for ho- 
mogeneous systems only, for which momentum is con- 
served, and elementary excitations s{k) for fixed wave- 
vector k are defined. It asserts that for drift velocities v 
such that \v\ < Vc = inf e{k)/k ,where k = \k \ , the flow 
is a superfluid, i.e., takes place without energy dissipa- 
tion[8]. For these systems the sound velocity Vs in the 
medium is defined by 
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A way to interpret and explain the Landau criterion 
is to view it as a spectral condition for a certain ap- 
proximate Hamiltonian. This spectral aspect has played 
recently an important role in the explanation of super- 
conductivity via superfluidity of the Cooper pairs in the 
BCS-BEC transition in traps [24]. At first one might be 
led to view the Landau condition as a (ground-state) sta- 
bility condition for the Galilei-transformed Hamiltonian 
Hij^N of Hn given by (2a): 



'Supported in part by CNPq. 



Hit_ 



N 



H, 



N 



Pn>0 



(2c) 



2 



(as long as \v\ < 77c), where P/v is the momentum op- 
erator: i.e., it is stabiUty (positivity) of the Hamilto- 
nian as viewed from the reference system attached to the 
pipe. Unfortunately, (2c) is not true ( see [22] and ref- 
erences given there) , which may be interpreted as a sort 
of metastability of the Landau state (see also [15]). (2c) 
is, however, true (for \v\ < Vc) when Hn is taken as an 
approximating Hamiltonian describing a finite number 
(0(1)) of noninteracting elementary excitations, of the 
type of those analysed in [10b] of the forthcoming model 
(21) (see also [15]), or the hamiltonian of the weakly in- 
teracting Bose gas (WIBG) introduced in ref.[9]: 

Hb = ^ Skaluk + 

X! ^["o«o(afeafc +a-fea-fe)] + 



where 



E^K«-feK)') + ii.c. 



2V 



{2d) 

Above, Ek denotes the kinetic energy, X/. is related to 
the Fourier transform of the interaction potential, V is 
the volume and a arc the usual Boson annihilation oper- 
ators, related to a cubic box with periodic boundary con- 
ditions, and the sums exclude the point k = [9]. When 
the Bogoliubov c- number substitution ao/VV = a.Q is 
performed on (2d), Hb becomes equivalent (after a Bo- 
goliubov transformation) to a Hamiltonian of (an infinite 
number of) independent elementary excitations, which 
satisfies (2c) [9]. It should be emphasized that (2d) may 
be controlled in several ways, both with and without the 
c- number substitution [9]. 

The explanation of the stability (2c) for a Hamiltonian 
describing independent elementary excitations is that it 
describes, of course, the idealized situation in scattering 
theory, in which the excitations are infinitely far apart 
and do not interact. Instabilities may be brought about 
by their interaction: this physical picture has been made 
rigorous for spin-waves in the infinite Heisenberg ferro- 
magnetic ground state [23] . It would be interesting to ex- 
tend the latter results to the Boson gas using the WIBG 
with the c-number substitution as Hamiltonian of the 
noninteracting quasi-particles. 

It is easy to see (see also [8]) that a nonzero critical 
velocity Vc implies 



Vs>0 
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and therefore a nonzero sound velocity is a necessary 

condition for the validity of Landau's criterion. This nec- 
essary condition has an analogue for inhomogeneous, i.e., 
non-translationally invariant systems such as trapped 
gases. Let the thermodynamic limit of the (ground state 
(g.s.), i.e., T = 0) compressibility kq be defined by 
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P = -{dE{N,L)/dV)N,T=o 



(7) 



is the g.s. pressure, and E{N,L) denotes the g.s. en- 
ergy of Htv- For homogeneous systems, and under quite 
general conditions on the pair interactions v it may be 
proved [6] that the limit in (5) exists and 
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and, by a macroscopic argument [7] , one expects that Vs , 
given by (4), is alternatively given by 
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On the other hand, (9) makes sense for general (in par- 
ticular inhomogeneous) systems such as (2). Comparing 
(9) with (4), we are led to adopt as necessary part of 
Landau's criterion for general systems the condi- 
tion: 
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In order to study (6) for the trapped gases (2) (with 
T ^ 0), we enclose the system in a large box and should 
consider the general limit on the r.h.s. of (6) where N, 
rather than Na/L, is fixed, so that wc arc not allowed to 
use the results for the GP limit (1), and the problem is 
difficult. One might argue intuitively that the trap po- 
tential Vt is "like a box", and enclosing the system in 
a large box whose volume is varied produces no change: 
thus P is zero in (7), and so is the inverse compressibility, 
yielding zero for the r.h.s. of (10). Thus the necessary 
part of Landau's criterion seems not to hold for trapped 
gases, but for a somewhat trivial reason. This is not so, 
however, because the trap potential is slowly varying, and 
very different from the steep ones which resemble a box. 
In the limit of very slow variation, one has the homoge- 
neous GP gas which has been thoroughly studied (see 
[12], chapter 5). The proofs of Bose Einstein condensa- 
tion (BEG) and superfluidity (according to the two-fluid 
criterion) in the latter reference have, in fact, been an ex- 
celent qualitative guide for the corresponding behavior of 
trapped gases, and for that reason we turn our attention 
now to the study of (10) in the homogeneous GP case. 

Let ^ — p denote the density. This may be fixed in the 
homogeneous GP limit but in a trap the particle density 
is inhomogeneous. We consider two versions of the GP 
limit which have been used: GPl: Na/L and 1/L fixed; 
GP2: p and Na/L fixed. 

Above, L and a are to be regarded as adimensional 
(i.e., divided by a length unit. In the trapped case, 
the natural length unit is the "trap extension" LqSc, see 
([12],p.47). As remarked in ([12], p. 40), GP2 may be im- 
plemented by replacing the interaction potential v in (2) 

by 
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with vi of unit scattering length, keeping vi fixed, with 
this choice a tends to zero as 



a = iV-2/3p-i/' 



since 
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and (1) holds. By scaling, the limit involved in the proof 
of ODLRO (off-diagonal long range order) [20] , which im- 
plies existence of BEC in a dilute limit, is equivalent 
to GPl ([12], Theorem 5.1, pg.40). The same happens 
with the proof of superfluidity by the two-fluid criterion 
(3) (see [12],pg.45). This equivalence does not extend, 
however, to the proof of (10), as wc shall sec. 
The g.s. energy in the thermodynamic limit. 
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is rigorously given in the dilute limit viz, 

pa^ < 1 , (15) 
,by the seminal result of Lieb and Yngvason [11]: 

eo(p) 



lim 

pa3^o Att fi p a 



1 



(6) 



The rigorous results referred to above (for a more com- 
plete account, see [12]) seem, so far, to support the as- 
sumption that the Bogoliubov expansion (in the forth- 
coming a = and (3 = 8(^ - v^)): 
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is an asymptotic series in region (15). We refer to 
([12],pg. 11) and Lieb's early review [10c] for details, 
references and history of (17), but remark that [13] is the 
simplest, albeit nonrigorous, approach to the derivation 
of the first two terms of (17). 

In the following considerations, the scattering length 
is a fixed number satisfying (17). In the limit (6) the 
derivatives with respect to V arc taken at fixed N and 
only thereafter the thermodynamic limit, hereafter just 

N 
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performed. We arrive ,by (6), (7) and (17), at 
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By (18) and (12) we obtain 
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Above, kJ^ — ^Tipp^/^ = ^p'^^'"^ denotes the compress- 
ibility of free fermions in three dimensions. We assume 
that f2{N), defined by (19b), is asymptotic of order two 

(in the parameter N~^) to the function g = TV^/^^^: 

this is a precise version of our assumption regarding 
asymptoticity. This yields by (19a) and (19b): 







(20) 



and thus the special form (10) of Landau's criterion is 
not satisfied for any regime of flow velocities. Thus GP2 
implies that the sound velocity is zero by (9). On the 
other hand, it is easy to see that (9) and (18) imply 
that GPl yields for the sound velocity the value 1/L 
(up to a constant factor with the dimension of velocity), 
i.e., a nonzero value! Both GPl and GP2 yield a value 
of the same order for the dilution parameter (15). 
Taking L = 1, GPl yields pa^ = N x (l/iV)^ = iV'^ 
(in the trapped case it is the mean density which en- 
ters (15), see [12]. pg.49) and by GP2 one finds also 
pa^ = p(7V-2/3)3 ^ o(iv-2). The values for Eq/N are, 
however, quite different: for GPl it is a nonzero constant 
proportional to Na, for GP2 it is , by (15) and (17), pro- 
portional to 7V~^/3 or V~^^^, which coincides with the 
behaviour of the free Bose gas with Dirichlet boundary 
conditions! This qualitative difference in behavior affects 
the leading term in (17), and therefore the compressibil- 
ity in (18) but, as we have seen, this does not matter for 
the presence or absence of superfluidity by the two-fluid 
criterion (3), since the latter holds even if Eq = 0. 

For homogeneous systems, wc view GP2 as the most 
natural, fundamental deflnition, and GPl as an equiv- 
alent condition for certain quantities. This is also indi- 
cated by the fact that the density p has no thermody- 
namic limit by GPl and is, in fact, the view adopted in 
[12], in which the GP limit is defined by GP2. From this 
point of view. Landau's condition should not be valid for 
the homogeneous GP case, and, we believe, also not for 
trapped gases in the GP limit, because, as argued above, 
the presence of the trap seems to strengthen the argu- 
ment in favor of a zero value for the compressibility. 

Condition (10) is, of course, equivalent to finiteness of 
the g.s. compressibility and is violated by the isother- 
mal compressibility of the free Bose gas below the tran- 
sition temperature (joining smoothly to T = 0). Clearly, 
the identity between the two expressions (4) and (9) is 
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highly nontrivial. To our knowledge, it has been rigor- 
ously proved (by E.H. Lieb [10b]) only for the model of 
N particles of mass m in one dimension with repulsive 
delta function interactions [10a] whose formal Hamilto- 
nian is given by 
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< < L, 1 = 1, 



(21) 



where c > 0, < i,j > denote nearest neighbours. We 
now explain why the compressibility of free fermions may 
be expected to appear in (19a). In the limit a — > im- 
plied by (12), (14) becomes a hard potential of zero range, 
just as the limit c ^ oo of (21) (Girardeau's model 
[14]). This seems to be the only model for which both 
the (double) spectrum of elementary excitations and the 
compressibility werc^ obtained rigorously without (as yet) 
unproved assumptions, the former in [10b]. The forth- 
coming derivation, which is not explicitly found in the 
literature but follows easily from [14], is elementary and 
provides a transparent physical reason for the validity of 
property (10). 

In the limit c ^ oo of (21), the boundary condition on 
the wave-functions reduces to 

ip{xi,. . . ,xn) =Q if Xj=xt 1 < j < £ < iV (22a) 

and the (Bose) eigenfunctions satisfying (22. a) simplify 
to 

(xi, . . . ,Xjv) = V'^ {xi, . . .,Xn)A{xi, . . .,Xn) 

(22b) 

where tp^ is the Fermi wave function for the free system 
of N particles confined to the region < Xi < L, i — 
1, . . . , N, with periodic boundary conditions: it automat- 
ically satisfies (22a) by the exclusion principle. In (22b), 

A{xi, . . . ,xn) = Y\_sgn{xj - xe) (22c) 

from the limit c ^ oo of [10a], or [14]. From (22b) and 
the fact that tpQ is nonnegative (the suffix zero referring 
to the ground state), it follows that Vo^ has constant sign 
in the A'^! regions in which configuration space is divided 
by the surfaces Xj = xe, and, thus, from (22. b), 



(22d) 



Since = 1 by (22. c), the correspondence (22. b) pre- 
serves all scalar products, and therefore the energy spec- 
trum of the system is the same as of the free Fermi gas, in 
particular the ground state energy equals (for N odd): 
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and V-'o' ^ Slater determinant of plane- 

wave functions labelled by wave-vectors 
ki,i= 1,. . . ,N, such that 
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with 
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the region (24a) being the "Fermi sphere" . A straightfor- 
ward analysis of the elementary excitations [10b] yields 
two branches, both corresponding to the same sound ve- 
locity by (4): 



Vs = 2K = np h/m 



(25) 



The r.h.s. of (9) follows by an elementary calculation: 
from the formula 

^ 1 

^ fc2 = - N{N + 1}{2N+1) 
fe=i 



we obtain from (8), and (22e): 
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By (26), Ko, given by (6), equals (the limit means A'' 
oo , with L ^ oo and N/ L = p): 



Ko 1 = lim ( 



(27) 



and thus the r.h.s. of (9) indeed equals the r.h.s. of 
(4), by (25). Thus in one dimension (27) shows that 

the ratio in (19b)equals one. Although the answer in 
higher dimensions is not known, it may be expected to 
be nonzero. Why does this ratio approach zero in the 
GP limit? We now attempt at an explanation, using the 
comparison with the homogeneous GP limit. 

It is the splitting in momentum space (18b), which is 
responsible for (10) in this model, with a if (oo) ^ 0, 
which itself is due to the (non-product) structure of the 
g.s. wave function (22b) - which would be identically zero 
if if (oo) = 0! In contrast, the latter condition is compat- 
ible with a g.s. wave function which is a product of plane 
waves in the same state (k = for periodic boundary 
conditions), as for the free Bose gas, leading to infinite 
compressibility. In the G.P. limit, the n-particle density 
matrices have a product structure ([12],pg.64), similarly 
to the free Bose gas, and quite in contrast to the rich 
structure of the correlation functions of the Girardeau 
model (see [21] for the case of Dirichlet and Neumann 
b.c). 

In ref. [15] we have formulated an alternate crite- 
rion applicable to rotational superfluidity, which depends 
on the property of ODLRO: the "macroscopic wave- 
fimction" associated to the latter property [20] is inde- 
pendent of the azimutal angle, for a system enclosed in a 
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rotating cylinder. The corresponding physical property 
is the "London rigidity" , which is shared by the free Bose 
gas [16]. In contrast, ODLRO is not required for Landau 
superfluidity, and, indeed, is proved to be absent for the 
Girardeau model in [17]. For dilute trapped Bose gases 
London rigidity was rigorously proved in [5], being thus 
one of the very few results in mathematical physics which 
have been verified experimentally ([18], Fig 7a, pg. 48). 

Although translational superfluidity of these gases was 
also rigorously proved in [5] using the two-fluid criterion, 
the necessity part of Landau's criterion (10) predicts oth- 
erwise, under assumptions on the asymptotic character 
of the Bogoliubov expansion. Thus question a) posed in 
the beginning seems to have a negative answer. Concern- 
ing question b), the issue can thus only be resolved by 
experiments, but the latter seem to be, as yet, inconclu- 
sive (see [19], Sect. 4.5, for a recent discussion of this 
point and references). In particular, if the compressibil- 
ity is very large (infinite by experimental standards), the 
two-fluid criterion is the relevant one in case the system 
is a (translational) superfluid, while the Landau criterion 
is the relevant one in case no translational superfluidity 
is found. 
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